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^ Abstract 

The possibility of an improvement of current NLO Monte Carlo event generators by 
means of choosing a suitable factorization scheme is studied. The optimal factorization 
t— I scheme for combining initial state parton showers and NLO hard scattering cross-sections 

is the ZERO factorization scheme, in which all NLO splitting functions vanish. However, 
pH it has turned out that the ZERO factorization scheme has a limited range of practical 

O | applicability. Hence, this paper is focused on searching for a factorization scheme which is 

applicable at the NLO in the full range of x relevant for QCD phenomenology and simul- 
taneously close to the ZERO factorization scheme (i.e. the corresponding NLO splitting 
functions are close to zero). 

1 Introduction 

^sO At present time many QCD cross-sections at parton level are known at the NLO accuracy and 

necessary algorithms, e.g. [H [21 13 HI [5] , for their incorporation in Monte Carlo event genera- 
tors have been developed. However, these algorithms attach initial state parton showers that 
are generated on the basis of the L0 splitting functions to NLO hard scattering cross-sections 
obtained in the standard MS factorization scheme because no satisfactory algorithm for gener- 
ating parton showers at the NLO accuracy has been found so far. The reasons why generating 
initial state parton showers at the NLO accuracy is difficult in the MS factorization scheme 
are following. The NLO splitting functions no longer correspond to basic QCD vertices, are 
expressed by much more complicated formulae than the L0 ones and are negative for some 
x, which prevents us from using straightforward probabilistic interpretation, which is crucial 
for Monte Carlo simulations. 

Initial state parton showers induce the scale dependence of parton distribution functions, 
which is described by the evolution equations. Monte Carlo simulations of initial state parton 
showers have some important advantages in comparison with the analytical calculations rep- 
resented by the evolution equations — they take into account transverse momenta, include 
the color description and the partons radiated in initial state parton showers participate in 
hadronization. 

Since current NLO Monte Carlo event generators attach L0 initial state parton showers 
to NLO hard scattering cross-sections corresponding to the MS factorization scheme, they 
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cannot predict some quantities, e.g. spectra of transverse momenta, at the NLO accuracyjj 
The deficiency lying in the combination of LO initial state parton showers and NLO hard 
scattering cross-sections calculated in the MS factorization scheme could be removed by using 
the ZERO factorization scheme, in which all NLO splitting functions vanish, and therefore 
all NLO corrections are included in hard scattering cross-sections. Attaching formally LO 
initial state parton showers to NLO hard scattering cross-sections is thus consistent if we use 
the ZERO factorization scheme. The change of the factorization scheme employed in NLO 
Monte Carlo event generators requires only to transform hard scattering cross-sections from 
the standard MS factorization scheme to a new one and to determine parton distribution 
functions in the new factorization scheme. The existing algorithms for parton showering and 
for attaching parton showers to NLO cross-sections need not be changed. 

The ZERO factorization scheme was studied in detail in [6j. Unfortunately, it has turned 
out that from the practical point of view, the ZERO factorization scheme cannot be applied 
at the NLO in the full range of x needed for QCD phenomenology. Hence, it is worth trying to 
find some factorization scheme which is applicable without any restrictions and simultaneously 
sufficiently close to the ZERO factorization scheme (i.e. the corresponding NLO splitting 
functions are sufficiently close to zero) . A factorization scheme satisfying these conditions will 
be called an approximately ZERO factorization scheme in the following. An approximately 
ZERO factorization scheme does not allow the construction of fully consistent NLO Monte 
Carlo event generators based on initial state parton showers that are taken formally at the 
LO. But its exploitation in NLO Monte Carlo event generators still makes sense because it 
can significantly reduce the deficiency caused by combining LO initial state parton showers 
with NLO hard scattering cross-sections calculated in the standard MS factorization scheme. 

This paper is focused on searching for an approximately ZERO factorization scheme. This 
search is described in detail in the next section. The summary and conclusion are presented 
in section [3j The notation used in this text is the same as that in [6]. 

2 Searching for approximately ZERO factorization scheme 

The NLO splitting functions should be as small as possible in an approximately ZERO fac- 
torization scheme. If NLO splitting functions are small, then they should have little influence 
on the evolution of parton distribution functions. Hence, the evolution of parton distribu- 
tion functions in any factorization scheme that can be considered as an approximately ZERO 
factorization scheme should be close to their evolution calculated with absolutely zero NLO 
splitting functions. This fact can be exploited in searching for an approximately ZERO fac- 
torization scheme. Since an approximately ZERO factorization scheme is aimed at using for 
the description of proton collisions, it is desirable to minimize the influence of NLO splitting 
functions on the evolution of parton distribution functions in the case of the proton. 

To quantify the influence of NLO splitting functions on the evolution of parton distribution 
functions in a factorization scheme FS, let us introduce a set of auxiliary parton distribution 
functions Do(x, M, FS, Mg) that are defined as follows: 

• their evolution in the factorization scale M is formally LO (i.e. the evolution is calculated 
with vanishing NLO splitting functions), independently of the factorization scheme FS, 

1 However, the NLO accuracy of some predictions, like those for total cross-sections, is not disturbed by 
using LO parton showers if the NLO parton distribution functions are used in the simulation of the hard 
scattering cross-section. 
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• the initial condition of the evolution is given as 
D (x, M = M s , FS, M s ) = D(x, M s , FS) 

where D(x,M, FS) represents the parton distribution functions of the proton. The auxiliary 
"zero" parton distribution functions Do(x, M, FS, Ms) are fully calculable from the parton 
distribution functions D(x, M, FS). In any factorization scheme FS that is close to the ZERO 
factorization scheme, the "zero" parton distribution functions Do(x, M, FS, Ms) should be 
close to the parton distribution functions D(x, M, FS) for all values of x and M, independently 
of the "starting" factorization scale Ms , which specifies the initial condition of the evolution 
of the "zero" parton distribution functions. 

An approximately ZERO factorization scheme was searched by minimizing the difference 
between B (x, M, FS, M s ) and D(x,M,FS). A small difference between D (x, M, FS, Ms) 
and D(x,M, FS) is a necessary, but not sufficient, condition for the factorization scheme FS 
to be close to the ZERO factorization scheme. The reason why this condition was used for 
searching for an approximately ZERO factorization scheme is that it takes into account the 
fact that an approximately ZERO factorization scheme is aimed at using for the description 
of proton collisions. 

In [6J, it has been shown that factorization schemes specified via the corresponding NLO 
splitting functions, which can be chosen at will, may have unexpected restrictions of their 
practical applicability — NLO splitting functions appearing at first sight as reasonable may 
specify a factorization scheme that cannot be applied in the full range of x relevant for QCD 
phenomenology. Unexpected restrictions of practical applicability are ruled out if the Mellin 
moments of the appropriate NLO splitting functions satisfy (see subsection 4.4 in [6]) 

4>) " P GG(n) ~ b) (P« (n) - P$(n,MS)) + 

+ P { Qkn) (P { Q ] Q {n) - („) + b) (P« („) - Pj»{n,m) - 
~ 2^S(«)^(«) (Pl&W - P G %) - P&(n,m + P^KMS)) = (1) 
for n for which 

b 2 -(P^(n)-P^)f-4P^(n)P§l(n)=0 and Ren > 1. (2) 

In the case of five massless quark flavours, the approximate values of n satisfying the preceding 
condition ^ are 1.9001 and 3.1798. The condition ^ depends only on singlet splitting 
functions, and therefore it does not put any constraints on the choice of non-singlet NLO 
splitting functions. Hence, the choice of non-singlet NLO splitting functions cannot cause any 
unexpected restrictions of practical applicability. Since an approximately ZERO factorization 
scheme is required to be applicable in the full range needed for QCD phenomenology, the 
influence of NLO splitting functions on the evolution of parton distribution functions should 
be minimized on a set of NLO splitting functions that satisfy the condition of practical 
applicability ([!]). 

The difference between Do(x, M, FS, Ms) and D(x,M, FS) was minimized on the set of 
factorization schemes in which the non-singlet NLO splitting functions vanish and the matrix 
of the singlet NLO splitting functions P^^(x) is expressed as 

TV 

p( 1 )(x) = p( 1) (x) + ^A fe Pi 1) (x) (3) 
fc=i 
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where the real parameters are completely arbitrary. All matrices in the preceding relation 
are square matrices of dimension 2. The set in which an approximately ZERO factorization 
scheme is sought is characterized by the functions (x), k > which have to be chosen 
such that the singlet NLO splitting functions P^\x) satisfy: 



• the condition of practical applicability Q, 




Both conditions must be satisfied independently of the values of the parameters The 
second condition, which is satisfied trivially in the ZERO factorization scheme, ensures that 
the sum of momentum fractions of all partons is independent of the factorization scale. The 
difference between Do(x, M, FS, Ms) and D(x,M, FS) was minimized for several different 
choices of the functions Pj^(x). 

The factorization scheme that has been found on the basis of the minimization of the 
difference between Do (a;, M, FS, Ms) and D(x,M, FS) for five massless quark flavours will 
be denoted as EPO in the following. The comparison of the parton distribution functions 
D(x,M, EPO) to the "zero" parton distribution functions Do (x, M, EPO, Ms) is displayed in 
figures [T] and [2j The non-singlet parton distribution functions are not displayed because 
any "zero" non-singlet parton distribution function q$ (x, M, EPO, Ms) is identical to the 
appropriate non-singlet parton distribution function q NS (x, M, EPO), which follows from the 
fact that the non-singlet NLO splitting functions vanish in the EPO factorization scheme. 
The parton distribution functions are compared at small factorization scales M because the 
differences between D(x,M, EPO) and Do (x, M, EPO, Ms) increase with increasing distance 
from the initial condition at the factorization scale MsJ^] The parton distribution functions 
in the EPO factorization scheme were calculated as follows. Using formula (4.5) in [2], the 
parton distribution functions were calculated at the factorization scale Mr = 10 7 GeV from 
the MS ones that corresponded to the MSTW 2008 set [7j. The parton distribution functions 
at a general factorization scale M were then obtained by the evolution from Mp to M in the 
EPO factorization scheme (the EPO factorization scheme is defined only for five massless quark 
flavours, and therefore the evolution corresponded to this number of massless quark flavours 
even for factorization scales smaller than the mass of the b quark) . Figures [T] and [2] show 
that the influence of the NLO splitting functions on the evolution of the parton distribution 
functions corresponding to the MSTW 2008 set [7] is almost negligible in the EPO factorization 
scheme. From the point of view of the evolution of the parton distribution functions of the 
proton, the EPO factorization scheme is thus very close to the ZERO factorization scheme. 

The comparison of the parton distribution functions D(x,M, MS) to the "zero" parton 
distribution functions Do(x, M, MS, Ms) is displayed^] in figures H and g The MS parton 
distribution functions that were used for the construction of these graphs corresponded to the 
MSTW 2008 set [7] for factorization scales greater than the mass of the b quark. The parton 

2 Since the evolution of initial state parton showers is backward in Monte Carlo event generators, the parton 
distribution functions are compared at factorization scales M smaller than Ms. 

^Similarly as in the case of the EPO factorization scheme, only the quark singlet and gluon distribution 
functions are displayed. However, contrary to the EPO factorization scheme, displaying the non-singlet par- 
ton distribution functions would make sense in the case of the MS factorization scheme because non-singlet 
parton distribution functions q NS (x, M, MS) differ from the appropriate "zero" non-singlet parton distribution 
functions qo S (x, M, MS, Ms). 
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Figure 1: Comparison of the quark singlet distribution functions in the EP0 factorization 
scheme. The differences between the curves are so small that the curves overlap. Note that 
the ranges of the axes of the left graphs differ from those of the right graphs. 
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Figure 2: Comparison of the gluon distribution functions in the EPO factorization scheme. 
Note that the ranges of the axes of the left graphs differ from those of the right graphs. 
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distribution functions at factorization scales less than the mass of the b quark were then 
obtained by the evolution from greater factorization scales in the MS factorization scheme for 
five massless quark flavours (and therefore the MS parton distribution functions differ from 
those of the MST W 2008 set [7j for factorization scales smaller than the mass of the b quark). 
Figures [3] and [4] show that the effect of the NLO splitting functions on the evolution of the 
parton distribution functions cannot be neglected in the MS factorization scheme, especially 
in the low x region. The influence of the NLO splitting functions on the evolution of the 
parton distribution functions also manifests itself by the dependence of the "zero" parton 
distribution functions Do (a;, M, FS, Ms) on the "starting" factorization scale Ms. Comparing 
figures [j}j4j one sees that the "zero" parton distribution functions Do (a;, M, FS, Ms) depend 
on Ms significantly more in the case of the MS factorization scheme than in the case of the 
EPO factorization scheme, for which the dependence is almost negligible. From the point of 
view of the evolution of the parton distribution functions of the proton, the EPO factorization 
scheme is significantly closer to the ZERO factorization scheme than the MS one. 

The singlet NLO splitting functions in the EPO factorization scheme are plotted in figures 
[5] and |6| The figures show that the singlet NLO splitting functions in the EPO factorization 
scheme are not small in comparison with those in the MS factorization scheme [8] and the 
singlet LO splitting functions. In the low x region, the singlet NLO splitting functions in 
the EPO factorization scheme even strongly dominate over the MS and LO ones. Hence, the 
EPO factorization scheme cannot be regarded as a factorization scheme close to the ZERO 
factorization scheme. The little influence of the singlet NLO splitting functions on the evo- 
lution of the parton distribution functions of the proton in the EPO factorization scheme is 
thus not a consequence of the closeness of the singlet NLO splitting functions to zero, but it 
is a consequence of their suitable shape. 

Besides the EPO factorization scheme, other factorization schemes in which the influence 
of the NLO splitting functions on the evolution of the parton distribution functions of the 
proton is almost negligible have been found. However, in all these factorization schemes the 
behaviour of NLO splitting functions and parton distribution functions is similar to that in 
the EPO factorization scheme. If the functions P^(a:) in formula |3) were chosen in such a 
way that they did not allow the singlet NLO splitting functions P^\x) to be large in absolute 
value in the low x region, then no NLO splitting functions with sufficiently small influence 
on the evolution of the parton distribution functions of the proton in the corresponding 
factorization scheme were found. It thus seems that in order that NLO splitting functions 
have little influence on the evolution of the parton distribution functions of the proton, which 
is a necessary condition that the factorization scheme specified by them be close to the ZERO 
factorization scheme, they have to be large in the low x region (provided that they satisfy the 
condition of practical applicability Q). But large values of NLO splitting functions are in 
contradiction with the closeness to the ZERO factorization scheme. The results of searching 
for an approximately ZERO factorization scheme based on the minimization of the difference 
between T)o(x, M, FS, Ms) and D(x, M, FS) thus indicate that there is no factorization scheme 
which has no restrictions of its practical applicability at the NLO and simultaneously is 
sufficiently close to the ZERO factorization scheme. This conclusion is also confirmed by the 
results of a direct minimization of the singlet NLO splitting functions ~P^\x) given by formula 
([3]) because the NLO splitting functions that have been obtained by the direct minimization 
do not have a negligible influcence on the evolution of the parton distribution functions of the 
proton, and therefore the corresponding factorization schemes cannot be considered as close 
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Figure 3: Comparison of the quark singlet and gluon distribution functions in the MS factor- 
ization scheme for large x. 
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Figure 4: Comparison of the quark singlet and gluon distribution functions in the MS factor- 
ization scheme. 
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Figure 5: Comparison of singlet LO and NLO splitting functions for five massless quark 
flavours. The singlet NLO splitting functions are plotted for the MS and EPO factorization 
scheme. 
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Figure 6: Comparison of singlet LO and NLO splitting functions at low x for five massless 
quark flavours. The singlet NLO splitting functions are plotted for the MS and EPO factor- 
ization scheme. Note that the singlet LO splitting functions and the singlet NLO splitting 
functions in the MS factorization scheme are multiplied by 100. 
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to the ZERO factorization scheme. 



3 Summary and conclusion 

A factorization scheme suitable for current NLO Monte Carlo event generators is the ZERO 
factorization scheme because its employment makes the combination of formally LO initial 
state parton showers and NLO hard scattering cross-sections consistent. But it has turned 
out [6] that the ZERO factorization scheme is not applicable in the full range of x needed for 
QCD phenomenology. The deficiency of current NLO Monte Carlo event generators lying in 
combining formally LO initial state parton showers and NLO hard scattering cross-sections 
can thus be removed by means of choosing a suitable factorization scheme only in a certain 
kinematic region in which the ZERO factorization scheme is applicable. Since this region 
covers, for instance, the production of heavy objects, which is important in searches for new 
physics, the use of the ZERO factorization scheme makes sense despite its limited range of 
practical applicability. However, if we require applicability in the full range relevant for QCD 
phenomenology, then the freedom in the choice of the factorization scheme does not allow us to 
remove the mentioned deficiency of current NLO Monte Carlo event generators because of the 
limited practical applicability of the ZERO factorization scheme. Hence, we have tried to find 
some factorization scheme which is applicable at the NLO in the full range needed for QCD 
phenomenology and simultaneously sufficiently close to the ZERO factorization scheme. Such 
a factorization scheme could at least significantly reduce the mentioned deficiency of current 
NLO Monte Carlo event generators. 

We have found NLO splitting functions which satisfy the condition of practical applicabil- 
ity and the influence of which on the evolution of the parton distribution functions of the 
proton is almost negligible. However, these NLO splitting functions are not small in compar- 
ison with those in the MS factorization scheme, and therefore the corresponding factorization 
scheme cannot be regarded as a factorization scheme close to the ZERO factorization scheme. 
Moreover, it has turned out that the condition of practical applicability ([I]) does not allow 
NLO splitting functions to be significantly smaller (in the absolute value) than those in the 
MS factorization scheme. The deficiency of current NLO Monte Carlo event generators lying 
in combining formally LO initial state parton showers and NLO hard scattering cross-sections 
can thus not be significantly reduced by means of choosing a suitable factorization scheme if 
we require applicability in the full range relevant for QCD phenomenology. 
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